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ABSTRACT 

With recent and archival Rossi X-Ray Timing Explorer (RXTE) X-ray measurements of the heavily 
obscured X-ray pulsar IGR J 16393-4643, we carried out a pulse timing analysis to determine the orbital 
parameters. Assuming a circular orbit, we phase-connected data spanning over 1.5 years. The most 
likely orbital solution has a projected semi-major axis of 43 ± 2 lt-s and an orbital period of 3.6875 ± 
0.0006 days. This implies a mass function of 6.5 ± 1.1 M Q and confirms that this INTEGRAL source 
is a High Mass X-ray Binary (HMXB) system. By including eccentricity in the orbital model, we find 
e < 0.25 at the 2er level. The 3.7 day orbital period and the previously known ~910 s pulse period 
place the system in the region of the Corbet diagram populated by supergiant wind accretors, and the 
low eccentricity is also consistent with this type of system. Finally, it should be noted that although 
the 3.7 day solution is the most likely one, we cannot completely rule out two other solutions with 
orbital periods of 50.2 and 8.1 days. 

Subject headings: X-rays: binaries — pulsars: individual (IGR J16393-4643) 



1. INTRODUCTION 

IGR J16393-4643 was initially discovered by Sugizaki 
et al. (2001) and is listed as AX J16390.4-4642 in the 
ASCA Faint Source Catalog. It was later detected by 
INTEGRAL during the first Galactic Plane Scan (GPS) 
of the Norma spiral arm (Bodaghee et al. 2006). The 
INTEGRAL GPSs have revealed many highly-absorbed 
sources that were not easily detectable with the soft 
bandpasses (< 10 keV) of most previous observatories. 
The brightest sources detected during the first year of 
the GPSs are listed in Bird et al. (2004). Among the 
sources, a few tens of them have never been detected 
prior to INTEGRALS observation. The second year cat- 
alog has doubled the number of sources that are either 
new or without a firm classification (Bird et al. 2006). 
Most of the new objects from the first catalog share com- 
mon characteristics such as their location in the Norma 
region and high intrinsic absorption (Lutovinov et al. 
2005). They are believed to be HMXBs (see Kuulkers 
2005), and IGR J16393-4643 [I = 338.0°, b = 0.1°) is 
probably a member of this class. That this system may 
be a HMXB has also been suggested by Sugizaki et al. 
(2001) and Combi et al. (2004) due to the huge hydro- 
gen column density towards the source, the hard spectral 
index (0.7-10 keV band), and its flux variability. 

Bodaghee et al. (2006) studied IGR J16393-4643 with 
INTEGRAL and XMM-Newton observations. XMM- 
Newton EPIC clearly detected pulsations at 2-10 keV 
with a 38 ± 5% pulse fraction and allowed for a refined 
position measurement; a potential counterpart, 2MASS 
J16390535-4642137, is about 2" away from the posi- 
tion. Using 15-40 keV INTEGRAL ISGRI data span- 
ning about 54 days, the pulse period was measured to be 
912.0 ± 0.1 s with a pulse fraction of 54 ± 24%. By fit- 
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ting a Comptonized emission model (comptt in XSPEC) 
to the non-simultaneous ISGRI and EPIC data, the col- 
umn density was measured to be Nn = (25 ± 2) x 10 22 
cm~ 2 . The spectrum also required an Fe Ka line, and 
marginally, an Fe K/3 line. The absorbed integrated flux 
was 4.4 x 10 -11 erg cm~ 2 s _1 in the 2-10 keV band and 
5.1 x 10 -11 erg cm -2 s _1 in the 20-60 keV band. By 
creating pulse phase-resolved spectra, Bodaghee et al. 
(2006) found that the pulsations affect the normaliza- 
tions but do not modify the spectral shape significantly. 

In this work, we use recent and archival RXTE ob- 
servations to determine the orbital parameters for IGR 
J16393-4643 using pulse arrival time analysis. We find 
three potential solutions to the orbit, although only one is 
particularly convincing. The observations are described 
in § 2. In § 3 we describe the pulse arrival time measure- 
ments and error analysis, and some basic features of the 
spectrum. The results of the pulse arrival time analysis 
are presented in § 4. We discuss the implications of our 
results in § 5. 

2. OBSERVATIONS 

Ten RXTE observations were obtained (Obs. ID 
91080) at four day intervals during 2005 October- 
November (hereafter; epoch 2), with each observation 
roughly 12 ks in duration (three to four RXTE orbits). 
Four additional 5 ks observations (two RXTE orbits) 
were obtained at 1 and 2 day intervals in 2006 February 
(epoch 3). In addition, we use an archival observation 
(Obs. ID 90069) spanning about 2.2 days from 1.2 years 
earlier (epoch 1). The epoch 1 data was contaminated by 
the recurring black hole transient 4U 1630-47, which is 
about 1° from IGR J16393-4643 on the sky. The pointing 
position for epochs 2 and 3 was offset by 0.2°, excluding 
4U 1630-47 from the field-of-view. A summary of these 
observations is shown in Table 1. Overall, these data 
extend over 1.56 yr and probe a range of timescales. 

Source and background light curves were created with 
data from the Proportional Counter Array (PCA; Jahoda 
et al. 1996) using standard FTOOLS. The PCA instru- 
ment consists of five identical multianode proportional 
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Fig. 1. — Barycenter-corrected and background-subtracted 3—24 
keV light curves. From top to bottom the plots are (a) a merged 
light curve of all observations, (6) epoch 1 data, (c) epoch 2 data, 
and (d) epoch 3 data. A characteristic error bar for each epoch 
is shown in the top-right corner of the panel. The sub-panels for 
epochs 2 and 3 show variability on shorter time scales. Note the 
range of time scales present and the higher flux level for the epoch 
1 observation due to contamination from 4U 1630-47. 



counter units (PCUs), operating in the 2-60 keV range, 
with an effective area of approximately 6500 cm 2 and a 
1° field-of-view at FWHM. PCUs and 2 were always op- 
erating, but any additional PCUs were also used if their 
good-time-intervals spanned entire RXTE orbits. Pho- 
ton energies were restricted to 3-24 keV, and the light 
curves were binned by 16 s, which balanced the desires 
for fine time resolution while maximizing the signal-to- 
noise. The SkyVLE 4 background model was used for 
the epoch 1 data, and the Faint background model was 
used for epochs 2 and 3. The arrival times of individual 
events were reduced to the solar system barycenter using 
the Jet Propulsion Laboratory DE-200 ephemeris (Stan- 
dish et al. 1992) and faxbary. The averaged PCU and 2 
light curves are shown in Figure Characteristic error 
bars for each epoch of data are shown in the top-right 
corner of the corresponding panel. 

3. PULSATIONS 

The pulse phase averaged flux varied by ~15% on hour 
time scales during epoch 1 (see Fig. [2>). Due to the 
contamination from 4U 1630-47, however, the variabil- 
ity cannot unambiguously be attributed to IGR J 16393- 

4 see http.heasarc.nasa.gov/docs/xte/pca_news.html 




Fig. 2.— Examples of IGR J16393-4643 pulsations using 16 s 
bins. The top two panels show what are considered to be clean 
pulsations, the middle two panels show rather noisy or irregular 
light curve segments, and the lower panels show examples of flares. 
The pulse arrival times for the six examples are superimposed as 
vertical dotted lines. Note how the maxima of the flares tend to 
occur near the maxima of the pulse profile. 



4643. During epoch 2, the phase averaged flux was con- 
stant at the 10% level for most of the observations, but 
varied by more than a factor of two over 3 hrs at one 
point. This outburst was in progress at the beginning of 
observation 91080-01-08 at HJD 53683.17. The epoch 3 
phase averaged flux was also constant within about 10% 
for the individual observations (two RXTE orbits), but 
varied by ~80% over the scale of days. For roughly half 
of the observations, the pulses are clearly defined, while 
at other times the light curves are noisy or are highly 
variable (flares) on ^100 s time scales. Figure |21 shows 
six examples of uninterrupted light curves, each spanning 
one orbit (hereafter called a "light segment" ) . 

3.1. Arrival Times 

Pulse arrival times were measured in a multi-step iter- 
ative process. We began by removing data within 500 s 
of a flare, which is defined to be when the flux is more 
than 75% above the mean level for the light segment. 
A preliminary pulse template was then created by man- 
ually aligning a selected set of epoch 2 light segments 
with clearly defined pulses. The alignment was accom- 
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TABLE 1 

RXTE Observation Log of IGR J16393-4643 



Obs ID 






Date 


rp- Q 

1 ime opan 


Hjxposure L ime 










(U.T.) 


(KSJ 


first 


^v^uuiiio/ b/ sr o u j 


90069-03-02 


1 


2004 


Jul. 


29.67-30.00 


28.52 


15.54 


41.49 b 




i 


2004 


Jul. 


30 00-31 99 


171.94 


50.22 


40.77 b 


Q1 flKD 01 01 

yiuou u j_-U-L 


2 


2005 


Oct. 


11.98-12.22 


20.88 


12.32 


10.89 


Q1 080-01-02 


2 


2005 


Oct. 


15.91-16.17 


22.08 


12.50 


11.50 


Q1 080 01 0'^ 


2 


2005 


Oct. 


20.05-20.28 


20.40 


11.63 


11.04 


91080-01-04 


2 


2005 


Oct. 


24.56-24.75 


16.86 


11.73 


11.47 


91080-01-05 


2 


2005 


Oct. 


28.44-28.66 


18.96 


12.19 


10.48 


91080-01-06 


2 


2005 


Nov. 


01.36-01.68 


27.48 


12.45 


10.20 


91080-01-07 


2 


2005 


Nov. 


05.30-05.65 


30.36 


11.50 


8.81 


91080-01-08 


2 


2005 


Nov. 


09.16-09.41 


21.30 


12.38 


14.70 


91080-01-09 


2 


2005 


Nov. 


13.03-13.28 


21.18 


12.35 


12.68 


91080-01-10 


2 


2005 


Nov. 


17.03-17.27 


21.12 


12.70 


10.25 


91080-01-11 


3 


2006 


Feb. 


13.80-13.90 


8.48 


4.69 


13.60 


91080-01-12 


3 


2006 


Feb. 


14.86-14.96 


8.63 


4.53 


10.41 


91080-01-13 


3 


2006 


Feb. 


16.83-16.95 


9.86 


5.09 


7.26 


91080-01-14 


3 


2006 


Feb. 


17.80-17.91 


9.26 


5.28 


11.68 



a For 3-24 keV photons. 

b Note the contamination during epoch 1 due to 4U 1630-47. 



plished by fitting a combination of sine functions to the 
selected light segments, and adjusting the phases accord- 
ingly. Manual alignment of the pulses is preferable to 
folding the entire light curve because it corrects for vari- 
ations in the arrival times due to the source's position in 
the orbit. All light segments were folded modulo 912.0 s 
(Bodaghee et al. 2006), and then the epoch 2 and 3 light 
segments were cross correlated with the preliminary pulse 
template. The cross correlation lag times were measured 
to the nearest second. The resultant pulse arrival times 
were then used to create a more refined pulse template 
including all epoch 2 and 3 data (minus the flares). We 
repeated this process until the pulse arrival times and the 
pulse template no longer changed. The final pulse tem- 
plate was then used to determine the pulse arrival times 
for the epoch 1 data separately. In total, sixty pulse 
arrival times were obtained. 5 They are presented with 
their associated errors in Table 2. The final pulse tem- 
plate, with a pulse fraction of 21 ± 1%, is shown in Fig. 
131 We note that the pulse template obtained from the 
PCA is fairly similar to the XMM-Newton folded light 
curve (Fig. 7, Bodaghee et al. 2006) with the main and 
secondary minima occurring at phases and ~0.3. 
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Fig. 3. — Final pulse template (3-24 keV) used to find the 
pulse arrival times via cross-correlation. The pulse fraction is 21 
± 1%. The dotted lines represent the phase corresponding to the 
pulse arrival times (Table 2). The error bars are lcr statistical 
uncertainties. 



3.2. Error Analysis 

Statistical errors on the pulse arrival times were calcu- 
lated using Monte Carlo simulations. For each light seg- 
ment, one hundred trial light curves were created. Each 
point on the light curve was chosen randomly from a 
Gaussian distribution centered on the measured number 
of counts, and with the standard deviation determined by 
Poisson counting statistics from source and background 
errors added in quadrature. Each trial light curve was 
cross correlated with the template, yielding slightly dif- 
ferent pulse arrival times. The quoted statistical uncer- 
tainty is the standard deviation of the difference in pulse 

5 We do not include two additional pulse arrival times from 
Bodaghee et al. (2006) due to ambiguity in the location of the 
pulse minima. 



arrival time relative to what was measured from each 
light segment. Statistical errors usually ranged from 3 to 
10 s, though the noisier light curves lead to errors of 15 
to 30 s, and a few were much larger. 

In addition to statistical errors, systematic uncertain- 
ties in the pulse arrival times can be caused by distor- 
tions in the average pulse profile. These distortions can 
be caused by varying flux, contamination from nearby 
sources (e.g., 4U 1630-47), and the energy dependence of 
the pulse shape can be affected by varying levels of local 
absorption. To estimate these errors, we aligned each 
folded light segment with the pulse template using the 
measured pulse arrival times, and each aligned and folded 
light segment was then used as a trial pulse template 
for all of the other folded light segments. The resultant 
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a variable overall multiplicative constant to adjust 
the normalization. The fourteen spectra were fitted 
simultaneously using XSPEC vll.3.2. Unfortunately, we 
found it difficult to constrain the spectral parameters, 
and especially the hydrogen column density, due to 
the Galactic ridge emission. From 3-5 (3-8) keV, the 
ridge emission accounts for about 75% (60%) of the 
total flux measured by the PCA. This explains the 21 
± 1% pulse fraction measured by the PCA, versus the 
38 ± 5% measured by XMM-Newton from 2-10 keV, 
and 54 ± 24% by INTEGRAL from 15-40 keV: The 
ridge emission effectively washes out the signal. Any 
uncertainty in the ridge emission normalization leads 
to even larger uncertainties in the inferred column 
densities. More importantly, since pulsar spectra are 
typically complex functions of both phase and energy 
(White et al. 1983), changes in the IGR J16393- 
4643 spectrum with pulse phase confuse the analysis. 

Table 2. Pulse Arrival Times 



Fig. 4. — Histogram of the difference in pulse minima versus the 
final template created by using all light segments pulse templates 
for all other folded light curves. The fitted Gaussian distribution 
has a = 5.76 s. 



arrival times from each trial generally created a Gaus- 
sian distribution centered on the measured arrival time, 
though occasionally there were outliers. Figure 0] shows 
a histogram of the difference in the measured pulse min- 
ima versus the final template, created by using all light 
segments as trial templates for all other light segments, 
for a total of 1770 total elements. A Gaussian was fitted 
to the histogram to find a systematic standard deviation 
of 5.76 s, and this error was added in quadrature to each 
statistical error. Although the statistical and systematic 
errors may be correlated, whereby addition in quadra- 
ture would overestimate the total errors, the histogram 
has statistically significant deviations from the normal 
distribution at a pulse minima difference of 20-30 s, sug- 
gesting that the errors associated with some of the pulse 
arrival times are underestimated in this scenario. More- 
over, if the pulse period is different from 912.0 s, folding 
each light segment on this period will result in slight 
smearing of the pulsations and will lead to additional 
systematic errors of order a few seconds. 

3.3. Spectral Analysis 

Motivated by the possibility of orbital modulation 
of spectral parameters, we extracted one spectrum 
from each epoch 2 and 3 observation. The epoch 1 
spectra were not analyzed due to contamination from 
4U 1630-47. The IGR J16393-4643 spectrum consists of 
emission from the pulsar, plus diffuse emission from the 
Galactic ridge due to the 1° field-of-view of the PCA. 
The pulsar emission was modeled with an absorbed 
power-law with a high energy cut-off plus a Gaussian to 
model iron line emission. The Galactic ridge emission 
was modeled using the Valinia & Marshall (1998) 
parameterization, consisting of an absorbed Raymond- 
Smith plasma component of temperature ~2-3 keV and 
a power-law component of photon index *~1.8. The 
Galactic ridge components were fixed at the Valinia 
& Marshall (1998) best-fit values, but we included 



Arrival Time 3, 


Statistical/Total 


Arrival Time a 


Statistical/Total 




Uncertainty 11 




Uncertainty b 


(s) 


(s) 


(s) 


(s) 





4.4/7.2 


39054559 


2.8/6.4 


5448 


8.2/10.0 


39377752 


3.1/6.5 


6372 


3.0/6.5 


39383220 


3.1/6.6 


11845 


72.0/72.3 


39395054 


4.9/7.6 


22749 


3.3/6.7 


39717343 


3.5/6.7 


56428 


7.6/9.6 


39734649 


3.8/6.9 


61892 


20.2/21.0 


39740112 


3.2/6.6 


79182 


5.0/7.6 


40068776 


5.9/8.2 


84670 


3.9/6.9 


40074244 


5.5/7.9 


90134 


3.6/6.8 


40079709 


3.6/6.8 


95582 


2.4/6.2 


40081529 


7.2/9.2 


101964 


2.8/6.4 


40391059 


1.4/5.9 


153000 


4.5/7.3 


40396526 


2.2/6.2 


169432 


14.2/15.3 


40401980 


4.3/7.2 


175807 


3.3/6.6 


40408369 


3.7/6.8 


182183 


3.5/6.7 


40725128 


2.8/6.4 


186610 


79.1/79.3 


40730603 


2.8/6.4 


192203 


15.2/16.2 


40736993 


2.3/6.2 


37956572 


4.4/7.3 


40742452 


3.9/7.0 


37962023 


2.6/6.3 


41070181 


3.5/6.7 


37973850 


3.5/6.8 


41075629 


5.0/7.7 


38296118 


3.0/6.5 


41081075 


2.5/6.3 


38302510 


5.2/7.7 


41087460 


4.5/7.3 


38307969 


3.6/6.8 


48740706 


2.2/6.2 


38314333 


1.8/6.0 


48746156 


4.0/7.0 


38658474 


3.0/6.5 


48832691 


27.5/28.1 


38664855 


3.3/6.7 


48838122 


355.1/355.2 


38670316 


2.5/6.3 


49002025 


25.9/26.5 


39043641 


3.5/6.7 


49006544 


23.8/24.4 


39049098 


1.8/6.0 


49086621 


4.3/7.2 



"Arrival times begin at HJD 53215.688912. 

b A systematic uncertainty of 5.76 s was added to each arrival time. 



A final source of confusion in the spectral analysis results 
from brief flaring episodes, which last for ^100-500 s. 
Interestingly, the strongest flares usually occurred at the 
pulse phase corresponding to the maximum of the pulse 
profile ((f) ~ 0.7 — 0.8, see Fig. |3J). To illustrate spectral 
distortion during the flares, we created a spectrum of 
the central 150 s of the flare seen in the lower-right panel 
of Fig. H(Obs. ID 91080-01-12-00), and compared it to 
the spectrum from the remaining portions of that RXTE 
orbit. The spectra were fitted simultaneously, with the 
hydrogen column density and Galactic ridge emission 
normalization tied to common values. For this particular 
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observation, a high energy cut-off was unnecessary, and 
7V H = 15+2 x i 22 cm -2 ( 90% connc i ence interval). The 
spectrum during the flare was best fit with a photon 
index T = 0.86 ± 0.14, while T = 1.46±j^a for the rest of 
the observation. This type of spectral hardening is not 
unexpected, yet it precludes constraining the hydrogen 
column density and power-law components simultane- 
ously due to correlations between these parameters. We 
leave a full phase resolved spectral analysis as future 
work. 



4. RESULTS 

The pulse arrival times were fitted with a seven param- 
eter model of the orbit and pulse period evolution. The 
arrival time of the nth pulse is given by 



t n — to + nP + —P P+a x sin i sin 



27r(t„ — 1 ) 

Pnrh 



where Pq is the pulse period at time to, P is the pulse 
period derivative and is assumed to be constant, a x sinz 
is the projected semi-major axis of the orbit, P Q rb is the 
orbital period of the system, and <j>o is the phase of the 
orbit at time to- The pulse number n is given by the 
nearest integer to 
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Fig. 5. — Lomb-Scargle periodograms for the epoch 1 data (top) 
and for the data from epochs 2 and 3 (bottom). Note the aliasing 
of the peaks in the lower panel due to the ~4 day separation of the 
epoch 2 observations. 



tn — ^0 _ t n — < 

(P) " P +0.5P(t n -t )' 



(2) 



We made sure the pulse period had not substantially 
changed from 912.0 s by creating Lomb-Scargle peri- 
odograms for the epoch 1 and epoch 2 /3 data separately, 
which are presented in Figure Determination of the 
orbital elements is complicated by the need to decouple 
the effects of the orbital Doppler delays and the intrinsic 
changes in the neutron star rotation rate. To minimize 
these concerns, we first searched for circular orbital so- 
lutions using the epoch 2 data only. The pulse arrival 
times were fitted to equation (1), using a fine grid of 
initial values for the pulse period, projected semi-major 
axis, orbital period, and phase. An extensive search led 
to five potential fits that can be seen as the minima in 
Fig. El which shows the lowest xt f° r circular orbits fits 
to the epoch 2 data as a function of pulse period. Given 
these fits, we then utilized all data by phase connecting 
the epoch 2 data to the data from epochs 1 and 3. To 
connect the three epochs of data, we used the best-fit 
circular orbit parameters for the epoch 2 data as initial 
values and cycled through a grid of pulse period deriva- 
tives [—3,3] x 10~ 8 s s^ 1 with a step-size of 5 x 10~ 10 s 
s _1 . One of the initial fits was rejected due to an inabil- 
ity to connect the data sets, and one was rejected due to 
a substantially larger value of xt, leaving three potential 
solutions. The results of the fits are presented in Table 
2] and the orbital pulse arrival delays are plotted in Fig- 
ures The error bars in the figures represent the 68% 
confidence level, and the secular change in the pulse pe- 
riods has been removed from the top and bottom panels 
so that the modulation is purely due to orbital motion. 

Of the three potential orbital solutions presented 
in Table ^ only orbital solution 1 appears convincing. 
This solution has a projected semi-major axis of 43 ± 2 
lt-s and an orbital period of 3.6875 ± 0.0006 days. The 
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Fig. 6. — Lowest Xv f° r an circular orbits to the epoch 2 data as 
a function of the pulse period (0.06 s bins). The pulse period was 
allowed to vary by 0.03 s in either direction, and its derivative was 
allowed to span the range P = [—3, 3] X 10 — 8 s s — 1 . 



X-ray mass function is 



fx(M) 



(M c 



(Mi + M c f 



4ir 2 (a x smi) 



6.5 ± 1.1 M 



0i 



(3) 

which is consistent with the HMXB interpretation pro- 
posed by Sugizaki et al. (2001), Combi et al. (2004), 
and Bodaghee et al. (2006). Finally, a pulse period of 
910.4 s with a derivative of —1.08 x 10~ 8 s s _1 implies 
a pulse period of about 911.3 ± 0.1 s at the time of the 
912.0 ± 0.1 s pulse period determination of Bodaghee et 
al. (2006). The 5cr difference potentially indicates that 
the pulse period derivative was not constant. This most 
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TABLE 3 

Circular Orbital Solutions to All Epochs 



Solution: 




1 


2 


3 


P i a 

1 pulse 


GO 


910.4178 ± 0.0006 


910.4028 ± 0.0007 


911.6167 ± 0.0005 


p 


(xl0~ 9 s s" 1 ) 


-10.81 ± 0.03 


-9.34 ± 0.04 


2.7 ± 0.3 


a x sin i 


(lt-s) 


43 ± 2 


60 ± 3 


238 ± 2 


^orb 


(d) 


3.6875 ± 0.0006 


50.2 ± 0.5 


8.1033 ± 0.0009 


MM) 


(M ) 


6.5 ± 1.1 


0.092 ± 0.014 


221 ± 6 


Xz/ 




2.55 


5.29 


7.72 


xl,rnm b 




1.00 


2.07 


3.02 



Note. — All errors are quoted at the 90% confidence level for a single parameter. The fits to eq. (1) each have 54 degrees 
of freedom. 
a Pulse periods at t = HJD 53654.98692. 

b The renormalized Xv was obtained by multiplying each total error by 1.6. 



likely was due to the average accretion rate increasing 
between 2003 February-March and 2004 July. 

Orbital solution 2 is statistically acceptable, but it is 
not convincing for four reasons: (1) The value of xl ( see 
Table ^) is twice as large as that of orbital solution 1. 
(2) The residuals to the phase connection to the epoch 1 
data systematically deviate from the orbital model. (3) 
There is a clear spread in the three to four pulse minima 
from each day of epoch 2 data, implying the pulse period 
and orbital period are not compatible. (4) Finally, with 
a mass function of 0.09 M Q , the system would have to 
be viewed nearly face-on (i < 13°) for the companion 
to be a high mass star (assuming a 1.4 M Q pulsar and 
a companion mass of 10 M©). We also note that about 
30% of the orbital phase is not sampled with these data. 

Orbital solution 3 is immediately suspect because its 
mass function is greater than 200 M Q , and a star of this 
mass is hard (if not impossible) to rectify with theo- 
retical predictions for the largest possible stellar mass 
(e.g., Krumholz 2005). If confirmed, solution 3 would 
represent an extremely important discovery because the 
orbital companion would more likely be an intermedi- 
ate mass black hole (IMBH) rather than a normal star. 
Since the pulsar is observable in X-rays there is obvi- 
ously a mass-donating companion, and therefore the sys- 
tem would necessarily be a triple system, which begs 
the question of how such a system could be dynamically 
formed. Moreover, there may be problems fitting the or- 
bit of the pulsar and mass-donating companion outside 
of the Roche lobe of the IMBH. Unusually large compan- 
ion mass aside, orbital solution 3 can be phase-connected 
over all three epochs, and a second orbit added to eq. (1) 
could possibly bring xt to a level comparable to the other 
solutions. On the other hand, the secular pulse period 
derivative is positive, which is clearly inconsistent with 
the 912.0 i 0.1 s pulse period measurement. Similar to 
orbital solution 2, nearly half of the orbital phase is not 
sampled with these data. 

Orbital solution 1 stands out as the most likely one, 
yet with xt = 2.55, the errors associated with the pulse 
arrival times would necessarily be underestimated by an 
average of about 60% to yield xt = 1-00- This is not 
unreasonable considering the quality of the data (see Fig. 

and the fact that there are statistically significant 
deviations from the Gaussian fit to the histogram used 



to estimate the systematic errors. As a final attempt 
to improve the fit, we added eccentricity to the orbital 
model by including the term 



— a x sin « sin 



Pnrh 



+ 4>o - u} p 



(4) 



in eq. (1), where e is the eccentricity and lo v is the lon- 
gitude of periastron. The additional term represents the 
first-order term in a Taylor series expansion in the eccen- 
tricity and is reasonably accurate for a mildly eccentric 
orbit (Levine et al. 2004). An F-test showed that the ad- 
dition of eccentricity to orbital models 2 and 3 was not 
significant, with chance probabilities of 0.16 and 0.09. 
Orbital solution 1, on the other hand, improved with a 
3.4% probability of the improvement occuring by chance 
(X 2 /^ |circ-*occ= 137.7/54 -> 120.9/52). With the im- 
proved fit to orbital solution 1, the best-fit parameters 
remained approximately constant with the exception of 
the projected semi-major axis which increased from 43±2 
to 55 ± 2 lt-s, yielding a new mass function of 13.1 ± 1.0 
Mq. The best-fit eccentricity was 0.15 ± 0.05. 

5. DISCUSSION 

By all indications, IGR J 16393-4643 is a HMXB sys- 
tem. It is heavily absorbed, with roughly an order of 
magnitude more absorption than can be attributed to 
the Galactic interstellar medium (Dickey & Lockman 
1990), indicating intrinsic absorption within the binary 
system. Such strong absorption is probably due to the 
stellar wind of a young companion star. Low-mass X- 
ray binaries, on the other hand, are generally not heav- 
ily extincted like IGR J16393-4643. In addition, sources 
like IGR J 16393-4643 that undergo occasional flaring are 
typically found in wind-fed systems. Still, the nature of 
the companion star for this source is not certain. The 
most likely orbital solution has an X-ray mass function 
of about 6.5 M Q , although the reader should bear in 
mind that including mild eccentricity to the orbit in- 
creases the mass function to 13.1 ±1.0 M Q . Nevertheless, 
it's position on the plot of pulse period versus orbital 
period — the so-called Corbet diagram (Corbet 1984) — 
implies that the companion is a supergiant with an un- 
derfilled Roche lobe rather than a Be star. The absence 
of transient behavior and the probable circular orbit fur- 
ther suggest that this is a supergiant system and not a 
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Be X-ray binary. Until recently, only 10 supergiant X- 
ray binaries were known to exist in the Galaxy (e.g., Liu 
et al. 2000), although many new members of this class 
are being uncovered by the INTEGRAL GPSs (Walter 
et al. 2006). Eclipses are often observed in pulsars with 
supergiant binary companions if the inclination of the 
system is high enough. By folding the PCA light curves 
modulo the orbital period of the three solutions, we do 
not find any evidence for eclipses, which allows the con- 
straint R c tan i < 18.5 R© if we assume the companion 
mass is much larger than that of the pulsar. Assuming a 
companion radius of 10 R©, the inclination must be less 
than about 60°. Nevertheless, further observations are 
required to obtain an unambiguous orbital solution for 
IGR J16393-4643. When combined with an optical mass 



function and a constraint on the inclination, the X-ray 
mass function would determine the masses of the two 
stars, which is important for constraining neutron star 
equations of state, and for understanding the evolution 
of binary systems. 

We would like to thank Riidiger Staubert for providing 
insight into the nature of the problem at the early stages 
of this work, and Jean Swank and Evan Smith for help 
in planning the RXTE observations. This research has 
made use of data obtained through the High Energy As- 
trophysics Science Archive Research Center Online Ser- 
vice, provided by the NASA/Goddard Space Flight Cen- 
ter. 
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Fig. 8. — Same as Fig. |7| but for orbital solution 2. 
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Fig. 9.— 



Same as Fig. Q but for orbital solution 3. 



